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ABSTRACT 

We propose that TeV 7-ray emission from blazars is produced by collisions near 
the line of sight of high energy jet protons with gas targets ("clouds") from the broad 
emission-line region (BLR). Intense TeV 7-ray flares (GRFs) are produced when BLR 
clouds cross the line of sight close to the black hole. The model reproduces the observed 
properties of the recently reported very short and intense TeV GRFs from the blazar 
Markarian 421. Hadronic production of TeV GRF from blazars implies that it is ac- 
companied by a simultaneous emission of high energy neutrinos, and of electrons and 
positrons with similar intensities, light curves and energy spectra. Cooling of these elec- 
trons and positrons by emission of synchrotron radiation and inverse Compton scattering 
produces delayed optical, X-ray and 7-ray flares. 

Subject headings: BL Lacertae objects: general, BL Lacertae objects: individual (Mrk 
421, Mrk 501), gamma rays: theory, radiation mechanisms: non-thermal 



1. INTRODUCTION 

The Fred Lawrence Whipple observatory has re- 
cently reported (Gaidos et al. 1996) the detection 
of two dramatic outbursts of TeV 7-rays from the 7- 
ray blazar Markarian 421 (Mrk 421). The first one, 
on May 7, 1996, had a doubling time of about one 
hour during which its flux increased above the quasi- 
quiescent value by more than a factor of 50. The 
second outburst, on May 15, 1996, lasted approxi- 
mately 30 minutes during which its flux increased by 
a factor of 20-25. These reports followed previous re- 
ported detections by the Whipple observatory of very 
strong bursts of TeV 7-rays from Mrk 421 which did 
not seem to be accompanied by a similar enhance- 
ment in GeV 7-ray emission (Kerrick et. al 1995; 
Macomb et al. 1995; Schubnell et al. 1996) and from 
Markarian 501 (Quinn et al. 1996). Although more 
than 50 active galactic nuclei (AGN) have been de- 
tected before in GeV 7-rays by the Energetic Gamma 
Ray Telescope Experiment (EGRET) on the Comp- 
ton Gamma Ray Observatory (see e.g., von Montigny 
et al. 1995; Thompson et al. 1995), all belonging 
to the blazar type of AGN. Of aU the EGRET 7- 
ray blazars, only the nearest, Mrk 421 at redshift 
z=0.031, has been detected in TeV 7-rays (Punch et 
al. 1992; Lin et al. 1992; Macomb et al. 1995; Gai- 
dos et al 1996). Also the blazar Mrk 501 at z=0.034, 
which is below the level of detectability by EGRET, 
has been detected in TeV 7-rays (Quinn et al 1996). 
It has been suggested, that perhaps all 7-ray blazars 
emit TeV 7-rays, but the opacity of the intergalactic 
space to TeV photons due to e+e~-pair production 
on the infrared background photons prevents us from 
seeing them in TeV photons (e.g., Stecker et al. 1993). 
Therefore, it is natural to assume that the two clos- 
est blazars, Mrk 421 and Mrk 501, are not unique and 
they well represent TeV 7-ray blazars. Although Mrk 
421 seems to flare in TeV 7-rays quite often, most 
of the EGRET blazars seem to show, within observa- 
tional limitations, less variability over very short time 
scales in the 30 MeV - 30 GeV energy range (see, how- 
ever, Mattox et al. 1997). 

The observed GeV and TeV 7-ray emissions from 
blazars are usually both interpreted as produced by 
inverse Compton scattering of highly relativistic elec- 
trons from their jet on soft photons, internal or exter- 
nal to the jet, (e.g., Maraschi et al. 1992, Bloom and 
Marscher 1993, Dermer and Schlickeiser 1993; 1994, 
Coppi et al. 1993: Sikora et al. 1994; Blandford 



and Levinson 1994; Inoue and Takahara 1996). Al- 
though the radio, X-ray and MeV-GeV 7-ray emis- 
sions are naturally explained by synchrotron radia- 
tion and inverse Compton scattering of high energy 
electrons in the jet, there are inherent difficulties in 
explaining TeV 7-ray emission as inverse Compton 
scattering of soft photons by highly relativistic elec- 
trons or positrons in pure leptonic jets. The main 
difficulty is fast cooling of electrons and positrons by 
inverse Compton scattering when they are accelerated 
to TeV energies in the very dense photon field of an 
AGN. Here we would like to propose an alternative 
model for TeV emission from blazars based on the as- 
sumption that AGN jets consist of normal hadronic 
matter (e.g., Mannheim and Biermann 1992). Tev 7- 
rays are produced efficiently by the interaction of the 
high energy protons in the jet with diffuse gas targets 
of sufficiently large column density that cross the jet. 
Such targets may be atmospheres of bloated stars, 
stellar winds, or gas clouds in the broad emission- 
line region (BLR) around the AGN. We show that 
the simple properties of hadronic production of high 
energy 7-rays, which are well known from lab ex- 
periments, together with the properties of the BLR 
of AGNs which are known from optical, ultraviolet 
and X-ray studies, explain both the observed quasi- 
quiescent emission and the outbursts of TeV 7-rays 
from blazars. We also predict the prompt emissions of 
TeV neutrinos, and delayed optical. X-ray and MeV- 
GeV 7-ray emissions that accompany TeV GRFs. 

2. THE HADRONIC COLLIDER MODEL 

We assume that 7-ray blazars are AGNs with 
highly relativistic jets of normal hadronic matter that 
point in the observer direction. Coulomb coupling of 
electrons to protons overcomes the Compton drag in 
its very dense photon field and makes it possible to ac- 
celerate the jet particles to very large Lorentz factors, 
r = l/\Jl - 13^ » 1. For r > 1, the kinetic energy 
of the jet resides mainly in protons. This energy is 
converted quite efficiently into TeV 7-rays which are 
beamed towards the observer hy pp ^ w^X; tt" 27 
when "clouds" with high column density from the 
BLR that surrounds the central region cross the jet 
near the line of sight. The quasi- quiescent emission 
is due to jet interactions with many, relatively dis- 
tant, gas "clouds" in the BLR. Strong GRFs are pro- 
duced when "clouds" cross the line of sight at much 
closer to the central engine. By "clouds" we mean 
diffuse material in the form of atmosphere of bloated 
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stars (Alexander & Netzcr 1994), stellar winds or real 
gas clouds. Hadronic production of TeV 7-rays is ac- 
companied by a simultaneous emission of TeV neutri- 
nos, electrons and positrons mainly via pp n^X; 
TT^ H^i^n', /Lt^ e^VeVn. The subsequent cool- 
ing of these electrons and positrons by synchrotron 
radiation, inverse Compton scattering and annihila- 
tion in flight, produces optical photons. X-rays and 
MeV-GeV 7-rays. 

3. HADRONIC PRODUCTION OF GAMMA 
RAYS 

The cross section for inclusive production of high 
energy 7-rays with a small transverse momentum, 
cpT = Et < 1 GeV in pp collisions (e.g., Neuhof- 
fer et al. 1971; Boggild and Ferbel 1972; Ferbel and 
Molzon 1974) is well represented by 



(l/27rpT)e-^-/^« fix), (1) 



where E « nipT is the incident proton energy, ain « 
35 mb is the pp total inelastic cross section at TeV 
energies, Eq ~ 0.16 GeV and f{x) ~ {l — x)^/y^ is a 
function only of the Feynman variable x = E^/E, and 
not of the separate values of the energies of the inci- 
dent proton and the produced 7-ray. The exponential 
dependence on Et beams the 7-ray production into 
< Et/E ~ 0.17/r along the incident proton di- 
rection. When integrated over transverse momentum 
the inclusive cross section becomes a^^da/dx « /(x). 
If the incident protons have a power-law energy spec- 
trum, dFp/dE w AE~", then, because of Feynman 
scaling, the produced 7-rays have the same power-law 
spectrum: 



"^^^ '^''-dE^NpCJir.IAE-", 



dE dEj 



(2) 

where is the column density of the target and / = 

Jlx'^-^f{x)dx. 

4. THE BROAD EMISSION-LINE REGION 

Detailed studies of broad optical and ultraviolet 
emission-lines, whose atomic physics is well under- 
stood, have been used to obtain detailed information 
on the BLR of AGNs. From their line-shapes, rela- 
tive strengths and their time-lag response to the vari- 
ations with time of the central continuum source, it 
was concluded that the BLR consists of high column 



density broad emission-line clouds (BLCs) that move 
with very large random velocities in the BLR: 

(i) The size of the BLR has been estimated from re- 
verberation mapping of both Seyfert 1 galaxies (e.g., 
Peterson 1993) and quasars (e.g. Maoz, 1997), with 
typical lag times between 10 days for Seyfert 1 galax- 
ies and 100 days for quasars, respectively. Typi- 
cally, Rblr ~ 3 X 10'^^L44^ cm, where L = L44 x 
10^^ erg is the luminosity of the AGN in ionizing 
radiation. 

(ii) The column density and mean density of the 
clouds were estimated from the ionizing flux of the 
central source and the relative line strengths from 
the partially ionized clouds. Very high densities and 
column densities were inferred. Typical values are, 



A/p ~ 10^3-24 and Up 



10^° cm 3, respec- 



tively. 

(iii) For spherical clouds of uniform density N^. = 
(4/3)ncrc. Consequently, the radii of BLCs are typi- 
cally , Tc = 10^^ ri2 cm, with ri2 ^ 0.1 — 100. 

(iv) The velocity distribution of the BLCs has been 
estimated from the profiles of the broad emission 
lines. Their full widths at half maximum indicate 
typical velocities of a few 10"^ km extending be- 
yond 10'' km at the base of the lines. Rever- 
beration mapping have clearly established that the 
velocities are not a radial flow (Maoz 1997). They 
seem to be consistent with the expected velocities of 
clouds orbiting massive black holes, Vc ~ ^jGMjR w 
1.15xlOV^8/i?i6 cTOs-1 , whereM = MsxlO^M© 
is the mass of the black hole and R = Rie x 10^^ cm 
is the distance from the black hole. 

(v) The covering factor, i.e., the fraction of the AGN 
sky covered by BLCs, was estimated from the ratio of 
Lya photons emitted by the BLCs to the H ionizing 
photons produced by the central continuum to be, 
Cblr ~ 0.1 . 

(vi) The total number of BLCs in the BLR was es- 
timated from the sizes of the BLR and BLCs and 
the covering factor. Assuming Cblr ^ 1, one finds 

NbLR = {4/3)CBLRRlLR/r?, . 

The UV spectrum of Mrk 421 does not show BLR 
and vL^{UOO A) - 1.3 x lO'^^h'^ erg s"^ (Kinney et 
al. 1991). The BLR may be swamped by beamed 
UV power-law emission from the jet. Since the broad 
lines equivalent width is > 30, weaker than in other 
AGNs, the isotropic ionizing continuum should be < 
lO^^erg and thus we estimate that Rblr ~ 3 x 
10^^ for Mrk 421. 
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5. QUASI-QUISCENT EMISSION AND GRFS 

The hadronic collider model predicts that the TeV 
7-ray emission from blazars fluctuates with time and 
shows spectral evolution, even if the jet properties do 
not vary with time on short time scales. The exact 
properties of individual flares depend on many un- 
known parameters of both the clouds (their geometry, 
density distribution, speed and trajectory relative to 
the jet and line of sight) and the jet (opening angle 
9 jet J exact orientation relative to the observer, parti- 
cle composition and differential energy spectrum of its 
high energy particles as function of distance from the 
jet axis and along the jet). The general properties of 
the quasi-quiescent emission and the flares, however, 
can be estimated using some simplifying assumptions: 

Consider a conical jet of particles from a source 
that is incident on a cloud at a distance R from the 
source. Let b and denote their impact parame- 
ter and angle relative to the jet axis. For the sake 
of simplicity, let us assume that the observer is lo- 
cated at inflnity on the jet axis. Most of the 7-rays 
seen by the observer must arrive from impact pa- 
rameters smaller than the critical impact parameter 
be ~ REa/E-y < Rdjet because of the exponential 
dependence of their production cross section (eq. 1) 
on Et- The number of clouds with b < be in the 
BLR is NblrEq /4:E'^. A quasi-quiescent background 
is formed by jet-cloud interactions only if this num- 
b er is la rge, i.e., if E^ < Eeru = \/NblrEo/2 w 
VCblrXRw/t 12) TeV. In that case the jet produces 
a quasi-quiescent 7-ray flux of 

^ ~ CsLRNpainlAE-". (3) 

and the BLR acts as a target with an effective column 
density of CblrNp, as long as Ej > Eo/9jet (below 
this energy the produced 7-rays are not beamed ef- 
fectively towards the observer). For E^ > Eerit the 
BLR emission is expected to fluctuate considerably. A 
flare with a large intensity contrast ratio (= maximal 
intensity/ quasi- quiescent intensity) is formed when a 
cloud crosses the line of sight at a relatively small 
distance. If the radius of the cloud is larger than the 
critical impact parameter, i.e., Ve > REq/E-^, then 
when the cloud blocks the line of sight, the 7-ray flux 
at photon energies E^ > {R/vejEo ~ 1.6i?i6/ri2 TeV 
flares up with a maximum intensity, 

HF 

—L^N,air.IAE--, (4) 



where Np is the average column density of the cloud. 
Thus the maximal intensity contrast of TeV GRFs 
compared with the quasi- quiescent emission is Np/CsLRNp w 
10 — 100. The duration of TeV emission in such flares 
is of the order of the time it takes the whole cloud to 
cross the line of sight, i.e., 

Tgrf ~ 2relve ^ 1.7 X 10Vi2i?if Af-'/' s . (5) 
The mean time between such strong flares is 

At « fcRFRBLRlCBLRhe « Q.f>L]J^^ C^lE^lyr^^T^ day 

(6) 

where Cblr = O.lCo.i and Tgrf = lO^Tg s. For 
E^ < 1.6i?i6/ri2 TeV the maximal GRF intensity 
is reduced by (reE^/REo)^ and the duration of the 
GRF is approximately the time it takes the cloud to 
cross the beaming cone: 

Tgrf ~ 2REo/veE^ ~ 3 x 10^ E^lyR^'^M^^^^ s. 

(7) 

Hence the GRF has the following general behavior 
when a cloud crosses the line of sight at a distance R: 
At energies below E^ ~ 1.6 x Ri(,/ri2 TeV, the in- 
tensity contrast increases with increasing energy while 
the duration becomes shorter. Above this energy both 
the intensity contrast and the duration become inde- 
pendent of energy. This behavior results in a spec- 
trum which becomes harder when the intensity in- 
creases and softens when the intensity decreases. The 
averaged quasi-quiescent emission spectrum therefore 
is softer than the spectrum of strong flares at peak 
intensity. 

The above predicted properties of the quasi- quiescent 

emission and the flaring of blazars in TeV 7-rays seem 
to explain quite well those observed for Mrk 421 and 
Mrk 501 (Punch et al. 1992; Lin et al. 1992; Kerrick 
et al. 1995; Macomb et al. 1995; Quinn et al. 1996; 

Gaidos et al 1996). 

6. NEUTRINOS FROM GAMMA-RAY BLAZARS 

Hadronic production of photons in diffuse targets is 
also accompanied by neutrino emission through pp 
TT^ — > /U^i'^ ; — > e'^UfjtVe- For a proton power-law 
spectrum, dFp/dE = AE~°' with a power index of 
a ~ 2, one flnds (e.g., Dar and Shaviv 1996) that the 
produced spectra of 7-rays and i/^'s satisfy 

dF^/dE RJ Q.ldF^/dE. (8) 

Consequently, we predict that 7-ray emission from 
blazars is accompanied by emission of high energy 
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neutrinos with similar fluxes, light curves and en- 
ergy spectra. The number of events from a 
GRF in an underwater /ice high-energy z^^ telescope is 
SNaTgrf J R^,ida,^,/dE^){dF,/dE)dE^dE,whcTcS 
is the surface area of the telescope, Na is Avo- 
gadro's number, <Ty^ is the inclusive cross section for 
Vf^p — > /iX, and i?^ is the range (in gm cm^'^) of 
muons with energy E'^ in water /ice. For a GRF with 
Fj ~ 10~^ cm~'^s~^ above = 1 TeV and a power 
index a — 2 that lasts 1 day, we predict 3 neutrino 
events in a 1 km^ telescope. Since the universe is 
transparent to neutrinos, they can be used to de- 
tect TcV GRFs from distant 7-ray blazars. If the 
reported GeV GRF from the brightest EGRET 7-ray 
blazar PKS 1622-297, which had a maximal flux of 
F-y ^ 1.7 X 10"^ cTO~2s~i photons above 100 MeV 
(Mattox et al 1997), was accompanied by a TeV GRF 
it could have produced 30 events within a day 
in a 1 km'^ neutrino telescope. 

7. X-RAY, MeV and GeV GRFS 

The production chain — !■ tt^ — > that 
follows jet-cloud collisions suddenly enriches the jet 
with high energy electrons. Due to Feynman scaling, 
their differential spectrum is proportional to the 7-ray 
spectrum 

due/dE « 0.5dn^/dE (9) 

and they have the same power-index a as that of the 
incident protons and the produced high energy pho- 
tons and neutrinos. Their cooling via synchrotron 
emission and inverse Compton scattering produces X- 
rays, and MeV and GeV 7-rays with a differential 
power-law spectrum 

dn-f/dE ~ (10) 

where (a -|- l)/2 « 1.5 — 2 . Hence, the emission 
of TeV 7-rays is accompanied by delayed emission 

of optical photons, X-rays, and MeV and GeV 7- 
rays. The peak emission of synchrotron radiation by 
electrons with a Lorentz factor Fg traversing a per- 
pendicular magnetic field H±(Gauss) moving with a 
Doppler factor d — {1 ~ (3cos6)/Th along the jet, oc- 
curs at photon energy (Rybicki and Lightman 1979) 
E't- - 5 X IO'^'^HaXIS keV. The electrons loose 
~ 50% of their initial energy by synchrotron radia- 
tion in 

Te « bxl(fV-^Hl^5 S « 1.2xl03jf-3/2^-l/2j-l/2 5. 

(11) 



Conseqiiently, the time-lag of X-rays is invc;rsely pro- 
portional to the square root of their energy. Similar 
time-lags for MeV-GeV 7-rays are expected if they are 
produced by Inverse Compton scattering from the self 
produced synchrotron photons. The integrated burst 
energy over the keV-GeV range is limited by the total 
electron energy to less than ^ 50% of the total energy 
in the TeV GRF. The spectral evolution of the X-ray 
flare (XRF) is a convolution of the spectral evolution 
of the high energy electrons and their cooling time. 
It is hardest around maximum intensity and softens 
towards both the beginning and the end of the flare. 
Because of electron cooling the spectrum should be 
harder during rise time than during decline of the 
flare. Indeed, all these features have been observed 
by ASCA (Takahashi et al. 1996) in the X-ray flare 
(XRF) that followed the TeV GRF from Mrk 421 on 
May 15, 1995. Detailed comparisons will be presented 
elsewhere (Dar and Laor 1997). 

8. SUMMARY AND CONCLUSIONS 

We have proposed an hadronic collider model to ex- 
plain TeV 7-ray emission from blazars. We have used 
simpliiying assumptions to derive from the model 
the main properties of quasi-quiescent emission and 
outbursts of TeV 7-rays, neutrinos and X-rays from 
blazars. The predictions agree with the Whipple, 
EGRET and ASCA observations of high energy 7- 
ray emission from Mrk 421 and Mrk 501. This seems 
to support an hadronic origin of TeV 7-rays emission 
from blazars. Although further optical, UV and X- 
ray studies of the broad emission-line region in AGN 
and observations of TeV 7 ray emission from blazars 
may provide more evidence for the hadronic nature of 
AGN jets, a decisive evidence will be provided by the 
detection of the predicted TeV neutrino fluxes from 
gamma ray Blazars. 
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